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Abstract 

At present, there are discrepancies between the measurements of several observables in 
B — > 7rK decays and the predictions of the standard model (the "B — > tcK puzzle"). 
Although the effect is not yet statistically significant - it is at the level of > 3a - it does 
hint at the presence of new physics. In this paper, we explore whether supersymmetry 
(SUSY) can explain the B — > ttK puzzle. In particular, we consider the SUSY model 
of Grossman, Neubert and Kagan (GNK). We find that it is extremely unlikely that 
GNK explains the B — > ttK data. We also find a similar conclusion in many other 
models of SUSY. And there are serious criticisms of the two SUSY models that do 
reproduce the B — > TtK data. If the B — > ttK puzzle remains, it could pose a problem 
for SUSY models. 
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Over the past several years, measurements have been made of a number of observ- 
ables in the decays of B mesons which are in disagreement with the predictions of the 
standard model (SM): e.g. indirect CP asymmetries in penguin-dominated B decays 
[1], triple-product correlations in B — > <pK* [2], polarizations in B — > ViV 2 decays (V* 
is a vector meson) [3], etc. None of these discrepancies is statistically significant, so 
that these disagreements only point to a hint of physics beyond the SM. Still, if these 
hints are taken together, the statistical significance increases. Furthermore, they are 
intriguing since they all point to new physics (NP) in b — > s transitions. 

Arguably, the most stringent discrepancy appears in B — > irK decays. Briefly, the 
effect goes as follows. There are four B — > ttK decays: B + — > tt + K° (designated as 
+0 below), B + -> n°K + (0+), B° d -> (-+) and -> 7r fT (00). In terms of 

diagrams [4], the amplitudes are given by 

A +0 = -P' , 

^2A 0+ = P' - T'e* - C'e* - P' EW , 
A -+ = p> _ T / e n ? 

^2 A 00 = -P' - P' EW - CV 7 . (1) 

In the above, we have neglected small diagrams and written the amplitudes in terms 
of the color-favored and color-suppressed tree amplitudes T' and C, the t-quark- 
dominated gluonic penguin amplitude P', and the color-favored electroweak penguin 
amplitude P' EW - (The primes on the amplitudes indicate b — > s transitions.) In addi- 
tion, we have explicitly written the weak-phase dependence (including the minus sign 
from V t lV ts [P']), while the diagrams contain strong phases. (The phase information 
in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix is conventionally 
parametrized in terms of the unitarity triangle, in which the interior (CP-violating) 
angles are known as a, (3 and 7 [5].) The amplitudes for the CP-conjugate processes 
can be obtained from the above by changing the sign of the weak phase (7) . Note that 
these diagrams include the magnitudes of their associated CKM matrix elements. 

The diagram P' EW is not independent. To a good approximation, it can be related 
to T' and C using flavor SU(3) symmetry [6]: 

P' EW = 3 - C -^R(T> + C>) + 3 - C -^^R(T> - C) . (2) 
4 ci + c 2 4 ci — c 2 

Here, the q are Wilson coefficients [7] and R = \(V t lV ts ) / (V* b V us )\. 

Now, in Ref. [4], the relative sizes of the diagrams were estimated to be roughly 

l:\P'\ , 0(A) :|n \ P ' EW \ , 0(X 2 ):\C'\. (3) 

where A ~ 0.2. With this estimate, the diagram C should also be neglected in the 
B — > ttK amplitudes above [Eq. (1)]. Note that the smallness of \C'\ is verified by 
more robust hadronic computations: \C'/T'\ ~ 0.3 is the prediction of NLO pQCD [8], 
and I C'/T'\ ~ 0.6 is the maximal SCET (QCDf) prediction [9, 10]. 
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There are nine measurements that have been made of B — > nK decays: the four 
branching ratios, the four direct CP asymmetries A l ^ p (ij = +0, 0+, — h, 00), and 
the mixing- induced CP asymmetry S°° P in B° d -> n°K° [11]. With this data and the 
expressions for the B — > 7T.K" amplitudes, one can perform a fit [12]. In the first fit, C 
was neglected in the B — > 7rlf amplitudes. A very poor fit was found: Xmin/ 'd.o.j '. = 
25.0/5 (1.4 x 10~ 4 ). (The number in parentheses indicates the quality of the fit, and 
depends on Xmin an d d.o.j . individually. It shows the percentage of the parameter space 
which has a worse Xmin- 50% or more is a very good fit; fits which are substantially 
less than 50% are poorer. 1.4 x 10~ 4 corresponds to a 3-4(7 discrepancy with the SM.) 
This result has led some authors to posit the existence of a " B — > nK puzzle" [13]. 

In the second fit, C was kept and the full amplitudes of Eq. (1) used. In this 
case, a good fit was found: Xmin/d-o.f. = 1.0/3 (80%). This has led some people to 
argue that there is in fact no B — > irK puzzle (for example, see Ref. [14]). However, 
\C'/T'\ = 1.6 ± 0.3 is required here. This is much larger then the theoretical estimates 
described above. If one takes this theoretical input seriously - as we do here - this 
shows explicitly that the B — > nK puzzle is still present, at > the 3cx level. 

The question now is: what type of new physics can explain the B — > nK puzzle? 
All NP operators in b — > sqq transitions take the form O l ^ q ~ sFibqFjq (q = u, d, s, c), 
where the F^j represent Lorentz structures, and color indices are suppressed. These 
operators contribute to the decay B — > nK through the matrix elements (nK\ O l ^' q \B). 
Each matrix element has its own NP weak and strong phase. Now, it has been argued 
that all NP strong phases are negligible [15]. In this case one can combine all NP 
matrix elements of B — > nK into a single NP amplitude, with a single weak phase: 

Y,(FK\O i >«\B)=A«j*« . (4) 

B — > nK decays involve only NP parameters related to the quarks u and d. These 
operators come in two classes, differing in their color structure: s a Fib a qpFjqp and 
SaFjbpqpFjqa (q — u,d). The matrix elements of these operators can be combined into 
single NP amplitudes, denoted A' ,q e t<s> 'i and A^^e^'i , respectively [16]. Here, $' and 
are the NP weak phases; the strong phases are zero. Each of these contributes 
differently to the various B — > nK decays. In general, A'' q ^ A' c ' q and $' g ^ $' g c . Note 

that, despite the "color-suppressed" index C, the matrix elements A' c ' q e are not 
necessarily smaller than the A ,q e l ®i. 

The B — > nK amplitudes can now be written in terms of the SM amplitudes to 
O(A) [P' EW and T' are related as in Eq. (2)], along with the NP matrix elements [16]: 

A +o = _p> + A >c,d e i*£ ; ( 5 ) 

v^A 0+ = P'-T' + P' EW + A' comb e^' - A c ' u e^ , 
A -+ = P > _ T , e n _ A>c,u e iK c } 

v^A 00 = -P' + P' EW + A' comb e^' + A c4 e^d , 

where A' comb e^' = -A> u e^ +A' d e i *'d. 
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In 1999, Grossman, Neubert and Kagan (GNK) proposed a new version of super- 
symmetry (SUSY) [17]. This model was promising for NP contributions to B — > ixK 
decays because it incorporates a new CP phase, and because it breaks isospin. In this 
paper we explore whether the GNK SUSY model can in fact explain the B — > nK 
puzzle, i.e. whether it gives the appropriate contributions to A.' ,comb e' 1 ®' , A^^e 1 ®'™ and 

We begin with a review of the GNK SUSY model, emphasizing those points which 
are important to our calculation. In R-parity-conserving SUSY models, the largest 
contributions to flavor-changing neutral current (FCNC) processes potentially come 
from the gluino-exchange SUSY box or penguin diagrams. The chargino and neutralino 
contributions are parametrically suppressed due to their small gauge couplings. The 
source of the gluino-mediated FCNC is the off-diagonal components in the scalar mass 
matrix in the basis where the quark mass matrices are diagonalized (super-CKM basis). 
Since we are interested only in the b — > s transition, we consider only the down-type 
scalar mass matrix. 

However, a generic form of scalar mass matrices is not acceptable because it leads to 
too-large contributions to FCNC processes (SUSY FCNC problem) and/or to the elec- 
tric dipole moments of the neutron and electron (SUSY CP problem). To evade these 
problems, people usually assume that SUSY is broken in a hidden sector and mediated 
to the observable sector by some flavor-blind interactions, such as gravity or gauge in- 
teractions. Then the squark mass matrices are diagonal matrices at a high-energy scale. 
The off-diagonal components in the squark mass matrices are generated by renormal- 
ization group (RG) running. In these popular models, such as minimal supergravity 
(mSUGRA) [18], anomaly-mediated SUSY breaking (AMSB) [19] or gauge-mediated 
SUSY breaking (GMSB) [20] models, the SUSY FCNC/CP problems are solved be- 
cause the RG-generated off-diagonal terms are typically very small and they do not 
include new sources of CP violation. On the other hand, as a consequence, they also 
cannot explain any possible deviation in the CP asymmetries in B decays. 

The GNK model assumes the following form of sdown mass-squared matrices: 



M~ 

d,LL(RR) 



( ™r (R)11 o o 



roj' , Mi =0 3X 3, (6) 



V m UR)v ^ 



L(R) 22 L(R) 2 ;i \ ' d,LR(RL) 

d,2 ,™jd,2 
L(fl) 32 " L HR) 33 



where off-diagonal components can be as large as the diagonal components. Although 
Eq. (6) is not supported by the above-mentioned popular SUSY-breaking models, it is 
well-motivated in SUSY GUT theories, where neutrinos are in the same supermultiplet 
as down quarks [21]. The zeroes in the above mass matrix are justified by the fact that 
the experimental results for K°-K° mixing, B®-B® mixing and B — > X s 7 are in good 
agreement with the SM predictions. In general they can get small non-zero values, but 
they do not affect our results much as long as we do not consider the very large tan f3 
region [22]. In our analysis below, we consider two scenarios: (i) only LL mixing is 
present (i.e. M~ i s diagonal), and (ii) both LL and RR mixing are present. 
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The mass matrix M- 



is diagonalized by 

d,LL ° J 



2 2 



i m 7 



(7) 



with 





/ i 










cos 6 L 




Vo 


— sin 9 L e 



-i8 L 




(8) 

Similarly, the exchange L <-> i? in (8) gives r fl . We restrict to — 7r/4 < < 7r/4 

= if RR mixing is absent) and — n < 5 L(R) < ir. 

The form given in Eq. (6) is not sufficient to give large SUSY contributions to P' EW . 
(Actually it is known that the gluino contribution to the Z-penguin is small [23].) 
In Ref. [17], the authors assumed that there is a significant mass splitting between 
the right-handed up and down squarks. Then the gluino box diagrams become the 
main source of the isospin breaking, and the scale a 2 s /m 2 SUSY (SUSY contribution) is 
comparable with a/M^ (SM contribution) 

We now turn to a review of the new-physics amplitudes A' ,comb e 1 ^' ; A' c,u e % ®'° anc [ 
A' c ' d e^d. These same NP amplitudes also contribute to B° s — > K + K and — > 
K°K°, and have been calculated within GNK SUSY in Ref. [24]. We closely follow 
this reference in our analysis, and use its treatment of the NP SUSY amplitudes. The 
color-allowed and color-suppressed NP amplitudes are given by 



A'' q e 1 ^ 
A'c,g e i%° 



x/2 

G_F 

V2 



-Xk(^ + c1)-(^ + cI) + (c! + ^1) 

-A,— c 8g (1 + — ) 



(9) 



where q — u,d. (Note: in Ref. [24], A' c ' q e l ®'° and A' ,q e l ®'i are switched.) In the above, 
*t = V£Vts and 



C?; Ci 



^1 



xM 
xM 



c 8g + 



c\ + 2cf 



= % m 



t(ml-ml)F B -(m 2 K )f K , 
1 2 ^ K 'ml)h, 



m- K )F» 



2ml 



m b (/2)(m q (i2) +m s (n)) ' 

2ml 

m b (/i)(m u (n) +m d (/j,)) 



(10) 
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where the c's and c's are Wilson coefficients of the effective operator in the GNK 
basis, m q is the averaged mass of up and down quarks, and naive factorization has 
been used for the hadronic matrix elements A vK and A K7I . Also, egg = —\tCs g , where 
C&g f — C 8g + C 5 in the standard basis. 

When only mixing between components 2 and 3 of the down-squark mixing matrices 
is allowed, the Wilson coefficients are given by 



r q 



a 2 sin 26 L e i&L \ 1 



F m~ 



18 



5 



a 2 sin 26 L e lt>L 



4V2G 



F m~ 



- (b L -> s L ) , 

7 1 3 2 

F { X br^ X 9 R g) + Q G ( X b L g' X VR~9) ~ Q^fe) ~ ^^Hg) 



(b L -> S L ) , 



alsm26 r e iSL \ 5 



Co = 



cl = 



AV2G F m 2 g 



(6l -> 5i) , 



a 2 sin 20^ fl ^ 3 

: ^ X b L gi X qL~9) + ( X bL9> X QL9) ~ 7, A \ X l 



4V2G F m 2 g 



0, 



3 v o L gi~iL,y/ ■ ^2 



(11) 

where x ab = ml/m 2 . Wilson coefficients with inverse chirality c's have exactly the 
same form, with the replacement L <-> R. Loop integrals are given by 



F(x,y) 
G(x,y) 
A(x) 
B(x) 



xlnx 



ylny 



(a: - - l) 2 (y -x)(y- l) 2 (x - l)(y - 1) ' 



a; 2 lnx 



y 2 Iny 



(x-y)(x-l) 2 (y-x)(y-l) 2 
1 (l + 2x)\nx 



(x-l)(y-l) 



2(1 -x) 6(1 -x) 2 ' 
ll-7x + 2x 2 lnx 



18(1 -x) 3 3(1 -a;) 4 ' 
Finally, for the chromomagnetic penguin, we have 



A, 



2a s e ff 8 a 2 sin (2# L )e 



iSr 



3tt 



where 



/f USY (*) = 



3 4 v / 2G F m| 

— 11 + 51x — 21x 2 — 19x 3 — 6x(l — 9x) log a; 
72(x- l) 4 



(12) 



(13) 



(14) 



The equations presented above allow one to calculate Wilson coefficients at the 
SUSY scale, taken to be m t . They then need to be renormalized to the scale \i = m^. 
The renormalization procedure described in Ref. [24] is used. This then gives the three 
NP SUSY amplitudes A'> comb e i *' , ^' 'V*«° and A' *?** at scale m b . 

We can now see if GNK can explain the B — > 7iK puzzle. In Ref. [12], fits were 
done with NP. The value of 7 was taken from independent measurements. (The value 
of 7 is the same as in the SM even in the presence of NP [25].) However, if all NP 
amplitudes are kept, there are more theoretical parameters (10) than measurements 
(9), and a fit cannot be done. For this reason, a single NP amplitude was assumed to 
dominate. Four possibilities were considered: (i) only A'' comb ^ 0, (ii) only A' c > u ^ 0, 
(iii) only A' c ' d ^ 0, (iv) A ,c ^e^'° = A' c ' d e^d , A'' comb = (isospin-conserving NP). 
A very good fit was found only if the NP is in the form of A'' comb e l ® (i.e. the SM 
electroweak-penguin amplitude). It is therefore often said that any NP invoked to 
explain the B — > ttK puzzle must contribute mainly to A' ,comb and little to A' c,u and 
A' c,d . (However, it should be noted that the fit with only A >c,u 7^ is not bad.) On the 
other hand, the GNK SUSY model gives nonzero values to all three NP amplitudes, 
and so the results of Ref. [12] do not hold. Another procedure must be used. 

Our analysis proceeds as follows. The three NP SUSY amplitudes depend on a 
number of theoretical inputs. We generate these randomly in the following ranges: 

• 300 < m-g < 2000 GeV, 

• 100 < rrig < 2000 GeV, 

• -tt/4 < 6 L , R < tt/4, 

• -7T < 5 LiR < 7T, 

• 7 = 67.631° [26], 

• m u ,m d (2 GeV) = 2.5 to 5.5 MeV [5], 

• m s (2 GeV) = 0.095 ± 0.025 GeV [5], 

• F B ^ K (q 2 = 0) = 0.34 ± 0.05 [27], 

• F B ^(q 2 = 0) = 0.28 ± 0.05 [27]. 

Note that we have taken = mg following SU(2) L symmetry. The weak phase 7 
is allowed to vary in the ±2<r range. For the other (theoretical) quantities for which 
an error is given, we take the range as ±l<r. With these values, A^ comb e % ^\ A' c,u e % ^' u 
and A' G ' d e l ^d are generated. 

Given the knowledge of the three NP amplitudes and 7, the B — > -nK amplitudes 
[Eq. (5)] and observables depend only on the two SM diagrams P' and T' (magnitudes 
and strong phases; P' EW is related to T'). We can therefore do a fit to see how well the 
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Mode 


RR[1(T 6 ] 


A-cp 


Scp 


B+ 




23.1 ± 1.0 


0.009 ±0.025 




B+ 




12.9 ±0.6 


0.050 ±0.025 




B° d 


-> TT-K+ 


19.4 ±0.6 


-0.097 ±0.012 




B° d 


^7T°K 


9.9 ±0.6 


-0.14 ±0.11 


0.38 ±0.19 



Table 1: Branching ratios, direct CP asymmetries A CP , and mixing-induced CP asym- 
metry S CP (if applicable) for the four B — > nK decay modes. The data is taken from 
Refs. [1] and [11]. 

B — > nK data is reproduced. If the Xmin * s acceptable, then we can conclude that the 
GNK SUSY model explains the B — > nK puzzle. If not, then it does not. 

In order to establish what constitutes an "acceptable" fit, we take our cue from 
ordinary observables. There, the 2a limit implies that 4.55% of the points of a Gaussian 
distribution lie outside this interval. In this spirit, we assume that the Xmin * s acceptable 
if the percentage of the parameter space which has a worse Xmin * s 4.55%, i.e. Xmin 
is taken to be < 11.31. (Note: in practice, there is no relation between Gaussian and 
Xmin distributions. We use the information from the Gaussian distribution only as a 
guide.) 

Before presenting the conclusions of this analysis, we must consider other con- 
straints. There are many constraints on SUSY models - electroweak precision tests, 
Am d , Am K , b — > s(d)~f, etc. However, by far the most stringent is that coming from 
B®-B® mixing. This is discussed in detail in Ref . [28] , and we closely follow the analysis 
presented here. Wefind that \Am s /Am s s M \ = 0.788 ± 0.195. This limits the SUSY 
contribution to B®-B® mixing. Using the expression given in Ref. [28], we compute the 
GNK SUSY contribution to |Am s |. To do so, three more theoretical parameters are 
needed, and we generate them randomly: 

• Bi = 0.86l8;g|[ [28], 

. B 4 = 1.1718$ [28], 

. B 5 = 1.94l°;i [28]. 

For each set of theoretical parameters generated, we check whether the constraint is 
satisfied (within ±2<r). 

The parameter space of GNK SUSY models is enormous - there are 12 SUSY 
parameters alone. In order to do our best to adequately sample this parameter space, 
500,000 sets of theoretical parameters were generated. For each set, we checked whether 
the B — > nK and the B®-B® mixing data were reproduced. The results are shown in 
Table 2, for the cases where (i) only LL mixing is allowed, and (ii) both LL and RR 
mixings are allowed. From this Table we see that the case with only LL mixing is 
preferred by the B° s -B° s mixing data. However, neither mixing scenario can explain the 
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Xlin < 11-31 


Am s 


both 


Xlin < 11-31 


Am s 


both 


74 


414357 


15 


102 


92844 


1 



Table 2: The number of points (out of 500000) which satisfy x 2 mi n( B ~* * K ) < n - 31 > 
the Am s constraint within ±2cr, and both constraints. In the left table, only LL mixing 
is allowed, while in the right table, both LL and RR mixings are allowed. 



B — > irK puzzle - in both cases, the B — > irK data is reproduced only in a tiny region 
of parameter space. The B®-B® mixing constraint reduces this (already small) region. 
We therefore conclude that it is very unlikely that the GNK SUSY model obeys the 
constraints from B — > ttK decays, and virtually impossible that it reproduces the data 
from both B — > ttK and B®-B® mixing. 

In Fig. 1, we present the SUSY contributions to several B — > nK observables. This 
helps identify which measurements lead to the large xlin f° r eacn °f the 500,000 GNK 
sets of parameters. In particular, we show R c vs. R n , where 



Rc = 2 
1 



R 



n 



2 



BR(B + -> tt°K + ) + BR(B~ -> ^K~) 
BR(B+ -»• n+K°) + BR(B~ -»• vr-i? ) 
"BR( J B° -> 7T-K+) + BR(5° -> 7T+K-) 
_ BR(B° d -> 7r°ir°) + BR(5° -> 7T°K ) 



(15) 



vs. S®? p , and vs. • All are scatter plots, showing the contribution of the 
GNK SUSY model to the various observables. As can be seen from this Figure, GNK 
has little difficulty in reproducing the combined R c and R n quantities. However, the 
SM can do this alone, showing that there is no discrepancy with the SM for R c and 
R n . GNK can also explain the A^? p and S®° p observables. Note that the SM alone 
has difficulty with these measurements. On the other hand, it is almost impossible for 
GNK to simultaneously reproduce A° CP and A~p~. This shows explicitly that it is the 
direct CP asymmetry measurements which are most problematic. 

There are several reasons that the GNK SUSY model cannot explain the B — * irK 
puzzle. First, for much of the parameter space, all three NP amplitudes are small. 
Thus, despite the presence of SUSY, the B — > itK system is basically described by 
the SM. However, we saw that the SM has a very poor fit in explaining the B — > nK 
observables, and so the same is true here. Second, the B — >• ttK measurements suggest 
that there is NP in the P' EW diagram (A' ,comb ). However, as indicated earlier, SUSY 
does not contribute significantly to P' EW . As a result, it is very difficult for SUSY to 
explain the B — > irK puzzle, and the GNK SUSY fits are generally poor. Third, we 
saw in the fits in which a single NP amplitude was assumed to dominate that the fit 
with A' c ' u 7^ was not bad. However, GNK generally does not generate only a large 
A' c,u - a large A'°' d is also usually found. Again, this leads to a poor fit. All of these 
can be seen in Fig. 2, which shows the plots of A' c,u vs. A l,comh and A' c,u vs. A' c,d . The 
bottom line is that it requires a very precise pattern of SUSY parameters to explain 
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Figure 1: Scatter plots of R c vs. R n (top), A^? p vs. S®° p (middle), and A^p vs. A~p~ 
(bottom), for LL mixing only (left), and LL and RR mixing (right). Horizontal and 
vertical lines represent experimental values within la and 2a. Plots include all 500,000 
GNK sets of parameters. Red points (dark grey in black and white) indicate only the 
SM piece of the SM + SUSY contribution. For A°® p vs. S®° p , the SM piece is a single 
dot because there is no direct CP violation when SUSY is not added. 
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Figure 2: Scatter plots of A' c > u vs. „4'> com6 (top) and A' c ' u vs. „4 /c ' d (bottom), for LL 
mixing only (left), and LL and RR mixing (right). Plots include all 500,000 GNK sets 
of parameters. 



the B — > TT-fT puzzle, and this is not found in most of the GNK SUSY parameter space. 

Of the very few points which satisfy both constraints, the great majority correspond 
to a large A' c ' u and a small A' c,d and A' ,comb . Also, all the points with Xmin(B ~^ < 
11.31 have a gluino mass less than 1.3 TeV. This is the only direct constraint on the 
SUSY parameters. 

As we have seen, it is extremely unlikely that the GNK SUSY model explains the 
B — > ttK puzzle. As noted earlier, there are other popular SUSY models: mSUGRA 
[18], AMSB [19], GMSB [20], etc, However they all automatically solve the SUSY 
FCNC/CP problems by not allowing any CP- violating phases. So these models cannot 
explain the B — > ttK data either. 

There are two SUSY models which do reproduce the B — > ttK data. They have (i) a 
large chargino contribution which allows large (2,3) mass terms in the up-squark sector 
[29] , or (ii) R-parity violation [30] . However, these two models have their own problems. 
The one with chargino contributions seems to be fine-tuned. It is not natural, i.e. it 
is hard to find a more microscopic theory which generates only (2,3) up-squark mass 
components in the LL or RR sector. And the R-parity-violating model lacks the beauty 
of SUSY, e.g. it does not have dark-matter candidates. We therefore conclude that if 
the B — > irK puzzle persists, SUSY models could have some difficulty 
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To summarize, the supersymmetry (SUSY) model of Grossman, Neubert and Kagan 
(GNK) [17] has great difficulty in explaining the B — > nK puzzle. The B — > ttK data 
can be reproduced in the GNK model, but only in a tiny region of parameter space. 
Other SUSY models, such as those with minimal supergravity [18], anomaly- mediated 
SUSY breaking [19] or gauge-mediated SUSY breaking [20], fare no better, as they 
do not allow any new CP-violating phases. There are two SUSY models which do 
reproduce the B — > ixK data [29, 30]. However, these models are either fine-tuned 
or lack some elements of ordinary SUSY theories. The B — > nK puzzle is still only 
a > 3<T effect, and so cannot be considered statistically significant. However, if this 
discrepancy with the SM remains in the years to come, it could pose a problem for 
SUSY models. 
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